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I. INTRODUCTION 



The infrared behavior of QCD Green functions is of fundamental interest. In addi- 
tion, these functions provide an important input for many calculations in hadron physics, 
for recent reviews see e.g. 0, 0). As infrared singularities are anticipated for some of 
these Green functions, non-perturbative continuum methods are needed to complement the 
knowledge gained in lattice Monte-Carlo calculations. Studies using different techniques, 
such as Dyson-Schwinger equations (see e.g. 0, 0, 11 and references therein), renormal- 
ization group methods Iff, stochastic quantization jjj, and lattice Monte-Carlo calculations 
(see e.g. (H H E3, 11, E| and references therein) have provided an unified picture of the 
infrared behavior of propagators in Landau gauge QCD in recent years. In this context the 
propagator of the Faddeev-Popov ghosts is of special interest: In the Landau gauge this 
propagator is infrared enhanced and diverges more strongly than 1/k 2 for k 2 — > 0. On the 
one hand, this reflects the Zwanziger-Gribov horizon condition 0, 0, On the other 
hand, in the Landau gauge it enforces the Kugo-Ojima confinement criterion [Til . 0. The 



accompanying infrared suppression of the gluon propagator relates to positivity violation for 
transverse gluons by imposing a cut in the gluon propagator 0. This resolves an old puz- 
zle already encountered in perturbation theory, which has led to the Oehme-Zimmermann 
superconvergence relations |l7| . 

These Landau-gauge studies are complemented by similar ones in the Coulomb gauge. 
Also in this gau ge t he infrared behavior of propagators is related to the Gribov problem 
and confinement |l8| . In addition, it has been shown that center vortices play a crucial role 
in the infrared enhancement of ghosts Thus, the following picture emerges: Degrees 
of freedom belonging to the indefinite-metric part of state space like the ghosts in Landau 
gauge or ghosts and Coulomb gluons in Coulomb gauge are infrared enhanced. This infrared 
enhancement is related directly to an effective cutoff at the first Gribov horizon 0, EH- 
The corresponding "excitations" are confining in that they mediate long-range correlations. 
Transverse gluons, on the other hand, are confined by these modes. The infrared part 
of the transverse gluon propagator is strongly suppressed. Besides an intuitive picture of 
confinement, this also provides a formal line of reasoning: Violations of positivity remove 
these states from the S matrix. Based on the relation of this picture to center vortices 
it seems natural to speculate about the importance of topological field configurations in this 
context. 

Although the picture, emerging from different methods described above, is in itself con- 
sistent and thus convincing, it is not yet complete. Finite-volume effects prevent lattice 
calculations to explore the extreme infrared. Functional continuum-based methods on the 
other hand necessarily involve truncations and the related errors are hard to control. For 
the functional methods the Landau gauge is advantageous due to its non-renormalization 



of the ghost-gluon vertex [2C| |21[ . To all orders in perturbation theory, the Landau gauge 
ghost-gluon vertex does not develop a genuine ultraviolet divergence, and especially for van- 
ishing incoming ghost momentum it stays bare. Furthermore, it has been argued that, in the 
extreme infrared, the gauge fixing term dominates over the Yang-Mills action |7| . Therefore, 
the infrared behavior of all Green functions is expected to be dominated by contributions 
involving ghosts. This hypothesis has been tested for the gluon and ghost propagators and 
has proven to be correct, thus alleviating very strongly the issue of truncation induced er- 
rors. At this point, a truly non-perturbative investigation of the ghost-gluon vertex has a 
twofold aim: First, it will add a further test of ghost dominance in the infrared. Second, and 
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more importantly, the result is crucial to assess the validity of recent investigations based on 



functional methods as all but the very first investigations 1 22J, |23j used a bare ghost-gluon 
vertex. Thus, we will present a semi-perturbative calculation of the ghost-gluon vertex based 
on its Dyson-Schwinger equation (DSE). For this project, it has proven advantageous that 
lattice results for the Landau gauge ghost-gluon vertex have been published very recently 
HI . We will compare our predictions to these data. 



This paper is organized as follows: To make it reasonably self-contained we will briefly 
discuss the non-perturbative gluon- and ghost propagators as they emerge from the solutions 
of their DSE's, truncated at the level of propagators. Then, a truncation for the DSE of 
the ghost-gluon vertex will be given. We then discuss the results of a semi-perturbative 
evaluation of this vertex. To this end, two types of input for the vertex to be calculated 
are used. This, and the comparison to lattice results, provides strong evidence that the full, 
non-perturbative ghost-gluon vertex is very close to the tree-level one for all momenta. 



II. GLUON- AND GHOST PROPAGATORS IN LANDAU GAUGE QCD 

Yang-Mills theory in D-dimensional Euclidean spacetime in the Landau gauge is described 
by the Lagrangian [2| 

£ = = \f; v f;„ + c a d,Dfc h (i) 

F% = d,Al - d„A; - g D r hc A\Al 



D a ; = S^ + gnf^A 



where F£ v denotes the field strength tensor, the covariant derivative, go the D 



dimensional gauge coupling, and / the structure constants of the gauge group. A® is 
the gluon field and c a and c a are the Faddeev-Popov ghost fields, describing part of the 
intermediate states of the gluon field. 

Within this framework, in Euclidean momentum space the Landau gauge gluon and ghost 
propagators, D fJill (p) and Da(p), can be generically written as 

D^(p,fi)=^6^--^ r j - 2 , D G {p,n) = - 2 , (2) 

where fi 2 denotes the renormalization scale, and Z(p 2 ,fi 2 ) and G(p 2 ,fi 2 ) are the gluon and 
ghost dressing functions. They can be determined from a solution of their DSE's P, H]| |26| 
using a well-established truncation scheme [i^, • A recent comparison of these solutions 
to the corresponding lattice results can be found in ref. ji^, H3- m the infrared, i.e. for 
infinitesimally small p 2 , these equations can be solved analytically 0, 0, and one finds 
simple power laws, 

z{p\ /i 2 ) - {p 2 r +2 - D ' 2 , g( P 2 , ~ (p 2 r« , (3) 

for the gluon- and ghost dressing function with exponents related to each other and to the 
dimensionality D. Here k is an irrational number, k ~ 0.595 for D = 4 and k ~ 0.398 for 



In these studies a ghost-gluon vertex which is an approximate solution to the corresponding Slavnov- Taylor 
identity has been employed. 
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FIG. 1: The truncated DSE for the ghost-gluon vertex. Dotted lines denote ghosts, and wiggly 
lines gluons. Lines with a dot indicate full propagators. Vertices with small black dots represent 
bare vertices and open circles represent full vertices. Contributions from the ghost-gluon scattering 
kernel have been neglected. 

D = 3 21, 2^1 . These analytical results depend slightly on the truncation scheme 2 |ill E^]. 
As mentioned above, this is in agreement with the Kugo-Ojima confinement criterion and 
Zwanziger's horizon condition. 

In four space-time dimensions the ghost- and gluon dressing functions can be used to 
define a non-perturbative running coupling ji^] 

a{p 2 ) = a^ 2 )G 2 {p 2 ^ 2 )Z{p 2 ^ 2 ) . (4) 

Due to the ultraviolet finiteness of the ghost-gluon vertex in Landau gauge, no vertex func- 
tion appears in this definition. Note that the r.h.s. of eq. Q is a renormalization group 
invariant, and thus a(p 2 ) does not depend on the renormalization point. From the analyt- 
ical power laws (J3J) one infers that the coupling has a fixed point in the infrared, given by 
cv(0) ~ 8.92/iV c . The infrared dominance of the ghosts imply that a(0) depends only weakly 
on the dressing of the ghost-gluon vertex and not at all on other vertex functions |2l|. 

In the following, for the calculation of the ghost-gluon vertex in four space-time dimen- 
sions the pointwise accurate fit [25| 



2 a(0) , 2 . cp 2K + dp iK 



,1k 



ln(e + aip 2a2 + b 1 p 2b2 ) ' ' 1 + cp 2K + dp 



Z{ * )= W)) Rip) ' Gip)= W)) ™- (5) 

will be used. It employs fitting parameters a 1; a 2 , 61,62 and c,d for the running coupling 
a(p 2 ) and the auxiliary function R(p 2 ), respectively. Here, 5 = —9/44, is the anomalous 
dimension of the ghost dressing function and a(fi 2 = (1.31 GeV) 2 ) = 0.9676. The six 
parameters of the fit are given by a x = 5.292 GeV _2a2 , a 2 = 2.324, 61 = 0.034 GeV~ 2b2 , 
b 2 = 3.169, c = 1.8934 GeV~ 2K and d = 4.6944 GeV~ 4K . For the calculation in D = 3, the 
numerical results for G{p 2 ) and Z(p 2 ) [26j are directly used. 



2 For D = 4 one can show, independent of any truncation, that k > 
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III. GHOST-GLUON VERTEX 



In the Landau gauge, the most general tensor structure of the ghost-gluon vertex with 
gluon momentum k and ghost momenta p and q is given by 

T abc (k; q,p) = ig D (q, (f abc + A abc (k 2 ; q 2 ,p 2 )) + k,B abc (k 2 ; q 2 ,p 2 )) , (6) 

where A abc and B abc are scalar functions describing the deviation from the tree-level form, 
and gn is the coupling constant. As there is no indication for a color structure differ- 
ent from the one occurring in perturbation theory [2] we assume that A abc (k 2 ; q 2 ,p 2 ) =: 
f abc A(k 2 ;q 2 ,p 2 ) and B abc (k 2 ; q 2 , p 2 ) =: f abc B(k 2 ;q 2 ,p 2 ). Note, that B is only relevant off- 
shell. 

The complete DSE for the ghost-gluon vertex is derived in the appendix 3 . According 
to the truncation scheme adopted for the propagators [lH |25|, the four-point function is 
neglected in the following. The truncated DSE is shown in Fig. ^ and is given in momentum 
space by 

T,(k;q,p) = rf\k;q,p) 

1 r dPoj 

xD G {u)r<®(q)D vX {u - q)T x (q -u;lu + k,p)D G (co + k) 

1 f dPoj 

-^9 2 D N C J j—^T^ p (k,uj,u + k) 

xD vX (u)rf\q)D G (u - q)T a {u + k;q- 0O,p)D p(T (u + k) , (7) 

where is the bare ghost-gluon vertex and r ^ vp the connected 3-gluon vertex. The 
momentum routing follows the same conventions as in ref. 0. 

Although the ghost-gluon scattering kernel is neglected, a self-consistent solution of this 
equation, together with the propagator equations, is of significant technical complexity. 
Fortunately, as we will see below, such a procedure is not necessary. It is sufficient to perform 
a semi-perturbative calculation, i.e. to do one iteration step in the ghost-gluon vertex DSE. 
If our starting hypothesis is correct, the resulting vertex should not significantly deviate from 
the input tree- level vertex. As a further test, we will also employ as input an ansatz for the 
ghost-gluon vertex, which is an approximate solution of the corresponding Slavnov- Taylor 
identity |22j |. 

r^;,,rt = ^(|g + |g-i). (8) 

First, we will present the results with the input vertices left bare and only the propagators 
dressed as described above. The results are displayed in Fig. |2 for the D = 4 case and in 
Fig. El for D = 3. Note that these functions have the proper ghost-antighost symmetry |2l| . 
e.g. A(k 2 ;q 2 ,p 2 ) = A(k 2 ;p 2 ,q 2 ) [3l[. The transverse part of the ghost-gluon vertex, 1 + A, 



is extracted employing 

P / k k 



q„ (<W - ^) r? c (£; q,p) = f abc (l + A(k 2 ; q 2 ,p 2 )) (9) 



The general structure of this DSE was already given in [3 
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FIG. 2: The normalized transverse part 1 + A (left panel) and the normalized longitudinal part 
B (right panel) of the ghost-gluon vertex for D = 4 in various kinematical regions. 
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FIG. 3: Same as Fig. |2 for D = 3. 



where A = q 2 k 2 — (q- k) 2 is a Gram determinant. The deviations of the transverse part from 
tree-level are clearly less than 20%. This is true for all momenta allowed by momentum 
conservation j3lj. In addition, also the longitudinal part, B {k 2 ; q 2 , p 2 ) , is smaller than 0.2 
for almost all momenta and finite everywhere. 

Thus, our results indicate that the full self-consistent solution will likely be very close to 
the tree-level form. A crucial further test is provided, if the non-trivial form (JHJ) is used as 
input on the r.h.s. of the the equation for the ghost-gluon vertex. 4 Several observations can 
be inferred from Fig. 0] First, also in this case deviations from tree-level are small, except 
for those kinematical regions where an infrared singularity is enforced by the ansatz (jSJ). 
Second, and even more important, the calculated vertex function A is much closer to the 
tree-level case than the input. A systematic study of possible input vertex choices yields the 
same result [3lJ. Furthermore, for D = 3 the results are very similar to the D = 4 ones |3l| . 

Finally, we want to compare our results to recent lattice results in Fig. El These 



4 The 3-gluon vertex is taken bare. Recent investigations showed that the self-consistently determined 
infrared divergence of the 3-gluon vertex does not change the result presented here [3^]. 
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FIG. 4: Same as Fig. [21 but here with eq. (JBJ) as input for the ghost-gluon vertex. The input 
vertices are denoted by (STI) and represented by crosses. 
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FIG. 5: The function 1 + A(0;p 2 ,p 2 ) for D=4 and two colors as compared to the corresponding 
lattice results The left panel shows the results obtained with symmetrically chosen momenta, 
the right panel with an asymmetric choice, see ref. [3| for further details. 



calculations have been performed for gauge group SU(2). Thus, we change the color prefactor 
of the loop diagrams in Fig. ^ accordingly. 5 Also, in the lattice calculation only the ghost- 
gluon vertex for vanishing gluon momentum has been determined, i.e. in our notation 
1 + A(0;p 2 ,p 2 ) has been calculated. Furthermore, the smallest momentum available on the 
lattice is 366 MeV, and this only at the expense of an asymmetrically chosen momentum, see 
ref. [24] for more details. Note that for symmetrically chosen momenta, see left panel of Fig. 
the lattice results for the vertex are less affected by the breaking of rotational symmetry. 
Given the systematic error in the lattice calculation, one can conclude that the lattice results 
are, within errors, consistent with the tree-level form at all momenta considered. Our results 
nicely match this behavior. In addition, we predict a slight decrease at small momenta. 

A feature of our results is the seemingly non-uniform limit for the functions A and B when 
all three momenta vanish. This is due to the ordering of the momenta when performing the 
limit. As these functions are finite, the full ghost-gluon vertex (including the corresponding 



5 This prefactor is f abc N c /2 for a general number of colors. 
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prefactors, see eq. ([fi])) is regular. In particular, for vanishing incoming ghost momentum, 
the ghost-gluon vertex is bare as expected (2(1 21 . 



IV. CONCLUSIONS 

We have presented approximate non-perturbative solutions for the ghost-gluon vertex in 
Landau gauge for Euclidean momenta in D = 4 and D = 3. To this end, we have employed 
non-perturbative gluon and ghost propagators. We used two types of input for the ghost- 
gluon vertex in the loop diagrams in order to estimate the behavior of a fully self-consistent 
solution. We have also compared our results to those of a recent lattice calculation. 

These results, when taken together, show rather conclusively that deviations of the ghost- 
gluon vertex from its tree-level value are very small, especially in the infrared. They thus 
validate the truncation scheme used to calculate the propagators of the Yang-Mills theory. 
More importantly, they confirm the strong evidence for infrared ghost dominance in Landau 
gauge, and thus for the Zwanziger-Gribov scenario, as they fulfill Zwanziger's hypothesis of 
a bare ghost-gluon vertex in the infrared. 

The results are not expected to change qualitatively when including matter fields, since 
the input propagators are rather insensitive 0, to quark contributions and no additional 
terms appear in the truncated DSE. This is quite distinct from similar calculations for 
the quark-gluon vertex, which find significant deviations from the tree-level form, possibly 
involving infrared divergences [34^ . 

In summary, the results presented here, nicely match a picture of confinement where the 
confining fields are on or near the Gribov horizon. They provide a further piece of evidence 
for a confinement mechanism of the Kugo-Ojima or Zwanziger-Gribov type. 
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APPENDIX: DERIVATION OF THE DSE FOR THE GHOST-GLUON VERTEX 

To derive the DSE for the ghost-gluon vertex it is only necessary to consider an action 
S g h that involves the contributions from ghosts in the Lagrangian ((I}, although due to 
the mutual coupling of the Green functions, the entire Lagrangian is implicitly relevant. 
Furthermore, we introduce J* a a and a a as sources for the fields A a , c a and c a , respectively, 
so that we can define the generating functional for full Green functions, 




(Al) 
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the generating functional for connected Green functions, W[J,a,a] = InZ, and the one for 
proper Green functions, T[A,c,c] = -W[J,a,a] + / d d x(A a ^ + a a c a + cV). The fields 
and sources are then given by 



SW 
5a a 
5T_ 

5c a 



= a 



5W 
5a a 
5T_ 

5c a 



= c 



6 -^ = A« 
ST 



5A 



- = T a 



(A2) 



where we use right derivatives for the following Grassmann fields, 

5 V 



5 5 
5c a '~ 5c a ' 



da a 6a a 

One way to approach the DSE for the ghost-gluon vertex is to start with the identity 

= J V[Acc]jJ^ exp (- J d D xC + j d D x (A a ^ + a a c a + c a a a )^j 
5S gh [A, c, c] 



(A3) 



5c b {y) 



+ <j\y) 



(A4) 



The expression in the brackets represents a full correlation function as generated by Z. 
Retaining non-zero sources, we now apply to the above expression the derivative 



5c c (z) 



5 2 Y 



5c d (v)5c c (z) 5a d (v) 



+ vanishing terms . 



(A5) 



Some terms vanish because the functionals can depend on pairs of Grassmann fields only. 
Throughout the calculation, great care is mandatory when dropping terms since most of 
them vanish only when setting sources to zero. We define the ghost and gluon propagators 
in position space 



ri=0 



5 2 W 
l Jx)fo b {y) 
5 2 W 



(A6) 
(A7) 



as well as the proper ghost-gluon vertex in position space 

5 3 T 



T abc (x; y, z) := — - 



5Al(x)8c b (y)5(f(z) 



(A8) 



r/=0 



After usage of the relation 
5(x-y)S° 



Sa a (x) 



J 5c d (z) 5a a (x) J ~ 5c d (z)5c b (y) 5a a (x)5a d (z) 



d D z- 



5 2 T 



5 2 W 



(A9) 
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one then obtains 
8 2 T 



8c b (y)8c c (z) 



Z[J,a,a] = d 2 8 bc 8(y - z)Z[J,a,a] 



+9Df b9h d y p j d D x 



8 2 Y 



8c d (v)8c c (z) 



A h p {y)c\y)Av)) • (A10) 



To find the DSE for the ghost-gluon vertex, we apply to equation (jA10|) the derivative 



8 



6A*{x) 



d D u 



5 2 r 



5A"(x)8Az(u) 5Jt(u) 



vanishing terms . 



(AH) 



One can now immediately set sources to zero to find the proper ghost-gluon vertex: 

8 3 Y 



g D f» h 8% J d D i 



8A a (x)8c d {v)8c c {z) 



[A h p (yy(y)^(v)) 



,1=0 



The decomposition of the full 4-point correlation function yields 

<^)^>)c%)c^))L f 



(A12) 



r?=0 



\ri=0 

5 2 W S 2 W 



5JWy)5JZ{u) 8a 9 {y)8a d {v) 
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8 3 W 



rv=0 



5 2 W 



5 2 W 



8JMy)8Jt{u) 8a 9 (y)8a d (v) 



77=0 



8 2 W 



8 2 W 



8JHy) 8Jt{u)8a 9 {y)8a d {v) 



8 3 W 



ri=0 



- I d n [rst) x 
8 3 T 



8Jz(u)8J*(r) 8a 9 (y)8a m (s) 8A\(r)8c m {s)8c n {t) 8J£(y)8a n (t)8a d (v) 

8 3 W 8HV 8 2 W 

' 8Jj}(y)8J^(u)8J^(r) 8a 9 (y)8a m (s) 8A\(r)8c m {s)8c n {t) 8a n (t)8a d (v) 



d D w 



8 2 W 



8 2 W 



8 2 W 



8Jt{u)8J k x (r) 8JHy)8Jl{w) 8a 9 (y)8 



a m {s) 



8 4 Y 



8 2 W 



x- 



8A l a (w)8A k x (r)8c m (s)8c n (t)8a n (t)8a d (v) 
8 2 W 8 3 W 8 3 T 



(A13) 



8 2 W 



8J^(u)8J k (r) 8JMy)8a 9 (y)8a m (s) 8A\{r)8c m {s)8c n {t) 8a n {t)8a d {v) 



T]=0 



The last term of the last line in eq. (jA13|) produces a 3PI-graph which, however, cancels 
in eq. (jA12|) with the first term. We now introduce two further definitions. The proper 
three-gluon vertex shall be denoted by 



(x, w, r) :- 



8 3 T 



8A a {x)8A^{w)8A{{r) 



(A14) 



ri=0 
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and the proper 4-point Green function involving two gluons and two ghosts is defined by 

5 A r 



TZ 9C (w,x;t,z):- 



5A^(w)5A a (x)5c 9 (t)Sc c (z) 



(A15) 
77=0 



After further decompositions of connected into proper 3-point correlation functions and 
using (jA9|) . eq. ()A12|) can be rewritten 6 as 

tt c (x;y,z) = g D f ahc dl8{y - x)5(y - z) 

+g D f hbm dl j d D [rstw]D h p 9 (y - t)D% n (y - w)f* nd (x; w, r)D d G e (r - s)Vf c (t; s, z) 

+g D f mbh dy J d D [rstw]D h G ?(y - t)D™(y - w)t^ d x (x, w, r)D d x l(r - ajf^a; t, z) 

-g D f mbh d 9 p j d D [tw}D h J(y - t)T™ 9C (w, x; t, z)D™(y - w) . (A16) 

The last step to take is to identify the bare ghost-gluon vertex which is derived from the 
Lagrangian (JTJ) as 

f ^ »• - S auA*m = 3Df ^ d " Av - x)Hy - z) ■ (A17) 

Thus one readily obtains 

g D f hbm d y p D h p i(y - t)D™(y - w) = J d d [uv}tf hbm ( U] y, v)D h p 9 (u - t)D™(v - w) . (A18) 

Using this, one can remove the spacetime derivatives in favor of bare ghost-gluon vertices 
and finally arrive at the complete DSE for the ghost-gluon vertex in position space: 

t abc (x;y,z) = f^ abc (x;y,z) 

+ J d d [rstuvw}D^ n (v - w)fl nd {x; w, r)D^ e (r - s)f f c (t; s, z)D h p 9 (u - t)f^ hbm (u; y, v) 

+ j d d [rstuvw]b™{u - w)f™ d x (x, w, r)D&(r - s)f ^(a; t, z)D h G 9 (v - t)T^ mbh (u; y, v) 

d d [tuvw]D 9 G h (v - t)f™ 9C (w, x; t, z)D™(u - w)tf mbh (u] y, v) . (A19) 

From this equation, eq. (JTJ) follows straightforwardly by Fourier transformation. 
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